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Abstract. 

Electroweak baryogenesis (EWBG) remains a theoretically attractive and 
experimentally testable scenario for explaining the cosmic baryon asymmetry. We 
review recent progress in computations of the baryon asymmetry within this framework 
and discuss their phenomenological consequences. We pay particular attention to 
methods for analyzing the electroweak phase transition and calculating CP-violating 
asymmetries, the development of Standard Model extensions that may provide the 
necessary ingredients for EWBG, and searches for corresponding signatures at the 
high energy, intensity, and cosmological frontiers. 
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1. Introduction 

Electroweak baryogenesis (EWBG) is one of the most attractive and promising ways 
to account for the observed baryon asymmetry of the Universe [H El |3]. As its name 
suggests, EWBG refers to any mechanism that produces an asymmetry in the density 
of baryons during the electroweak phase transition. While many specific realizations of 
EWBG have been proposed, they all have many features in common, and in this review 
we attempt to describe them in a unified way||] 

The initial conditions assumed for EWBG are a hot, radiation-dominated early 
Universe containing zero net baryon charge in which the full SU {2)lxU (l)y electroweak 
symmetry is manifest [HI |12l [131 [131 HI] • As the Universe cools to temperatures below 
the electroweak scale, T < 100 GeV, the Higgs field settles into a vacuum state that 
spontaneously breaks the electroweak symmetry down to its U{l)em subgroup. It is 
during this phase transition that EWBG takes place. 

Successful EWBG requires a first-order electroweak phase transition. Such a 
transition proceeds when bubbles of the broken phase nucleate within the surrounding 
plasma in the symmetric phase. We illustrate this process in Fig. [1} These bubbles 
expand, collide, and coalesce until only the broken phase remains. 

Baryon creation in EWBG takes place in the vicinity of the expanding bubble walls. 
The process can be divided into three steps [4]: 

1. Particles in the plasma scatter with the bubble walls. If the underlying theory 
contains CP violation, this scattering can generate CP (and C) asymmetries in 
particle number densities in front of the bubble wall. 

2. These asymmetries diffuse into the symmetric phase ahead of the bubble wall, where 
they bias electroweak sphaleron transitions [151 HB] to produce more baryons than 
antibaryons. 

3. Some of the net baryon charge created outside the bubble wall is swept up by 
the expanding wall into the broken phase. In this phase, the rate of sphaleron 
transitions is strongly suppressed, and can be small enough to avoid washing out 
the baryons created in first two steps. 

We illustrate these three steps in Fig. |2| 

These EWBG steps satisfy explicitly the three Sakharov conditions for baryon 
creation [IT]. First, departure from thermodynamic equilibrium is induced by the 
passage of the rapidly-expanding bubble walls through the cosmological plasma. Second, 
violation of baryon number comes from the rapid sphaleron transitions in the symmetric 
phase. And third, both C- and CP-violating scattering processes are needed at the phase 
boundaries to create the particle number asymmetries that bias the sphalerons to create 
more baryons than antibaryons. 

All the ingredients required for EWBG are contained in the SM. Unfortunately, 
EWBG is unable to explain the observed baryon asymmetry within the SM alone. The 

t See Refs. [1 [51 H EHi il [10] for previous reviews of EWBG. 
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Figure 1. Expanding bubbles of electroweak- broken phase within the surrounding 
plasma in the electroweak-symmetric phase. 

first impediment is that the SM electroweak phase transition is first-order only if the 
mass of the Higgs boson lies below rrih ^ 70 GeV [HI [19]. This is much less than 
the current experimental lower bound of rrih > 115.5 GeV [201 EI]- Even if the phase 
transition were first order, the CP violation induced by the CKM phase does not appear 
to be sufficient to generate large enough chiral asymmetries [221 [231 121] • 

Therefore an essential feature of all viable realizations of EWBG is new physics 
beyond the Standard Model (SM). This beyond the SM (BSM) physics must couple 
to the SM with at least a moderate strength, and it must be abundant in the thermal 
plasma at the time of the electroweak phase transition. Together, these two conditions 
imply the existence of new particles with masses not too far above the electroweak 
scale and direct couplings to the SM. Thus, a generic prediction of EWBG is that new 
phenomena should be discovered in upcoming collider and precision experiments. It is 
this property that sets EWBG apart from many other mechanisms of baryon creation. 

Because of the prospects for experimental probes of EWBG, it is particularly 
important to achieve the most robust theoretical predictions for the baryon asymmetry 
within this framework as well as for the associated phenomenological implications 
within specific BSM scenarios. Consequently, we review both progress in developing 
the theoretical machinery used for computations of the baryon asymmetry as well 
as developments on the phenomenological front. The former entail a mix of non- 
perturbative Monte Carlo studies and various perturbative approximations. Work on 
the phenomenological side includes applications to specific BSM scenarios, such as the 
Minimal Supersymmetric Standard Model (MSSM), and the delineation of consequences 
for collider studies, low-energy probes of CP-violation, and astrophysical observations. 

The plan for this review is as follows. In Section |2] we discuss the electroweak phase 
transition in greater depth, concentrating on its strength and other characteristics. Next, 
in Section |3] we describe in more detail the creation of asymmetries in the CP and baryon 
charges during the phase transition. Some of the ways the new ingredients required for 
EWBG can be studied in the laboratory are studied in Section |4j Finally, Section |5] is 
reserved for our conclusions. 
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Figure 2. Baryon production in front of the bubble walls. 



2. The electroweak phase transition 

Baryon creation in EWBG is closely tied to the dynamics of the electroweak phase 
transition (EWPT). In this transition, the thermal plasma goes from a symmetric state 
in which the full SU{2)l x U{1)y gauge invariance is manifest to a broken one where 
only the U{l)em electroweak subgroup remains [HI [T31 E]. As discussed above, the 
transition must be first-order and proceed through the nucleation of bubbles of the 
broken phase. In this section we will discuss the dynamics of this phase transition and 
describe the role it plays in EWBG. 

The transition from symmetric to broken phase in the SM can be characterized by 
the vacuum expectation value (VEV) of the Higgs field H = {H~^, H^)^ that transforms 
as (1, 2, 1/2) under SU{?>)c x SU{2)l x U{1)y- A field basis can always be chosen such 
that only the real component of develops a non-zero expectation value. Thus, we 
will write 



The symmetric phase corresponds to = and the broken phase to 7^ 0. Note that 
(in unitary gauge) the masses of the and Z° weak vector bosons and the fermions 
are proportional to 0. 

The features of this transition that are most relevant for EWBG are (a) its character 
(first order, second order, cross over); (b) the critical temperature Tc and the bubble 
nucleation temperature T„ that describe when it occurs; (c) the sphaleron transition rate 
Fsph that governs the rate of baryon number generation and washout; and (d) the bubble 
nucleation rate. These features have been studied using a broad range of theoretic tools. 

The most robust computations of many of these quantities are performed using 
non-perturbative, Monte Carlo methods. However, given the level of effort required 
to perform such studies, they have only been applied to a few specific theories of 




(1) 
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electroweak symmetry breaking. Instead, perturbative methods have been used much 
more frequently to study the dynamics of the EWPT in a broad range of BSM scenarios. 
Perturbative analyses can also provide helpful insight into some aspects of phase 
transition dynamics that be may be less accessible with Monte Carlo methods. It 
should be emphasized, however, that the application of perturbation theory to EWPT 
physics is fraught with uncertainties as well as the potential for ambiguities. In the SM, 
for example, one often finds the transition temperature computed perturbatively to be 
significantly lower than the value obtained from Monte Carlo studies for a given value 
of the Higgs boson mass. Nevertheless, given the widespread use of perturbation theory, 
we begin by reviewing this approach, first laying out the conventional treatment and 
then commenting on various difficulties. We subsequently review salient features of the 
non-perturbative analyses. 



2.1. Perturbative analyses 

In a perturbative analysis of the electroweak phase transition, the central object is the 
(renormalized) finite-temperature effective potential. This quantity coincides with the 
free energy of the cosmological plasma [H |25l [26l |27], provided it is reasonably close 
to thermodynamic equilibrium. The key feature of the effective potential is that the 
expectation value (p oi the Higgs field is that which minimises its value. 
To one-loop order, the effective potential is given by [8] 

K//(0,T) = K)(0) + yi(0) + Ayf)(0,T) , (2) 

where Vq = — m^0^/2-|- A0^/4 is the tree-level potential, Vi(0)''°'' is the one-loop effective 
potential at T = 0, and T) contains the leading thermal corrections. 

The expression for Vi(0) is well-known 



2si 



(3) 



where the sum i runs over all particles in the theory, each with rii degrees of 
freedom, field-dependent mass mi{(f)), and spin. Furthermore, we have assumed a mass- 
independent renormalization with /i as the renormalization scale and the Cj are scheme- 
dependent constants. Choosing fi ~ max{mi(0)} optimizes the perturbative expansion. 
Let us also mention that Fadeev-Popov ghosts are massless and decouple if we work in 
the Landau gauge = 0), in which case = 3 should be used for each vector boson 
and contributions from Goldstone bosons included as well. 

The thermal corrections are given at one-loop order by jS] 

AVr\^,T)^ J: n,I^^jM)- J: n,I^^jM) (4) 

i=boson ' j=fermion ^ 

where and Jj are loop functions. For small arguments, x <^ 1, they have the 
expansions [29] 

J,(a;^) = - ^ + - - ln(^V%) + 0{x^) (5) 
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with ln(a6) ~ 5.4076 and ln(a/) ~ 2.6351. At large arguments, a; ^ 1, both loop 
functions reduce to [29] 



This form shows the familiar Boltzmann suppression of particles much heavier than the 
temperature. 

To illustrate the effect of thermal corrections on the Higgs potential, it is helpful 
to write the potential in a simplified approximate form using the high-temperature 
expansions of Eqs.(5][6), noting that heavy particles (m ^ T) decouple quickly. This 
yields [29] 

T) ^ DiT' - T',)ct? - ET<P' + , (8) 

where D and A are slowly- varying functions of T (but not (f)). 

In the limit of = in Eq.(|8| the phase transition is second order, with a transition 
temperature of T = Tq and a Higgs expectation value for T < To of 



To\r-f{i-Tyn). (9) 



For non-zero E in Eq.(|8]), the phase transition becomes first-order. Starting from 
T ^ To, a second minimum away from the origin develops when T = Ti with 



where the temperature-dependent coefficients D and A are to be evaluated at T = Ti. 
This second symmetry-breaking minimum becomes degenerate with the origin at the 
critical temperature T^., and becomes deeper at lower temperature, as illustrated in 
Figj3j The degree to which the phase transition is first order is typically characterized 
by (pc/Tci where 0c is the location of the minimum at T^.. In terms of the parameters in 
the potential, it is 

^ = H. (11) 

Some time after the temperature falls below Tc, regions of the cosmological plasma tunnel 
to the deeper broken minimum and the phase transition proceeds by the nucleation of 
bubbles. 

Before discussing the dynamics of such a first-order phase transition, it is worth 
examining the validity of the perturbative expansion outlined above. This expansion 
is known to break down at very high temperatures when the thermal loop expansion 
parameter (yf^T^/m^(0) becomes large [HI 120] (where is the coupling entering in 
the loop). Indeed, we saw above that the leading thermal corrections, which are only 
generated by loops in this formalism, completely change the vacuum structure of the 
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Figure 3. Schematic temperature dependence of the effective potential 



theory and induce a restoration of symmetry at very high temperatures. The breakdown 
of the perturbative expansion can be postponed by resumming the most dangerous 
thermal corrections by incorporating thermal mass corrections in the propagators. The 
net result of such a daisy resummation is to generate an additional term in the effective 
potential [31j : 

Vt'-'^ = ~Y.^^ Hi^^T) - mm''' , (12) 
{h}' 

where the sum runs only over scalars and longitudinal vectors, and nf' is the field- 
dependent thermal squared mass: 

frf{<P) = m\<P) + Yi{T) , (13) 

with n(T) oc the thermal contribution to the mass. 

The daisy correction is particularly important for a first-order transition because it 
affects primarily the crucial cubic term. For example, suppose the contribution to the 
cubic term comes from a scalar with a zero-temperature mass of m^(</)) = gcfP' with a 
thermal correction of n(T) = kT^. The would-be cubic term becomes 

^E^' = ^9'''<P' ^ ^ [9<P' + ^T'Y^" (14) 

When n(r) is large relative to m^(0), this corrected expression ceases to behave as a 
cubic in and the phase transition might no longer be first-order. 

When the electroweak phase transition is first-order, it proceeds by the nucleation of 
bubbles of the broken phase within the surrounding plasma of symmetric phase. Bubble 
nucleation is governed by thermal tunnelling [32] from the local minimum at = to 
a deeper minimum at 7^ 0. In nucleating a bubble there is a competition between 
the decrease in free energy, proportional to bubble volume, with the increase due to the 
tension of the wall, proportional to bubble area. As such, there is a minimum radius for 
which a bubble can grow after it is formed, and this limits the tunnelling rate. Bubble 
formation and growth only begins in earnest when this rate exceeds the Hubble rate. 
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which occurs at some temperature T„ < Tc, called the nudeation temperature. Once a 
sufficiently large bubble is formed, it expands until it collides with other bubbles and 
the Universe is filled with the broken phase. The typical profile and expansion rate of a 
bubble wall can be computed from the effective potential [291 [331 [M] , taking into account 
frictional effects from scattering with surrounding particles in the plasma [33 ESI EZ] • 
As we will see in the next section, baryon creation processes are very sensitive to the 
speed and profile of the walls. 

A first-order electroweak phase transition is not sufficient for successful EWBG. 
The transition must also be strongly first-order. Within the context of perturbation 
theory calculations, the quantitative condition for a strongly first-order phase transition 
has typically been taken to be 

1^ > 1 ■ (15) 

This ratio approximates a factor that appears in the rate for sphaleron transitions in the 
broken phase within the bubble walls. When this condition is not met, these transitions 
will wash out the baryons created by EWBG. As we will describe below, the necessity 
of a very strong phase transition is one of the reasons why EWBG does not work in the 
SM, and is therefore a strong motivator of new physics beyond the SM. 



The condition of Eq.(15) is a frequently-applied approximation, but a precision 
calculation of the baryon abundance from EWBG requires a more careful analysis. For 
one, the relevant temperature for the dynamics is the slightly lower bubble nucleation 
temperature rather than the critical temperature, although it is often the case that 
Tn — Tc [Ml [39] • A more serious problem is the lack of gauge invariance. In particular, 
it is well-known that the vacuum expectation value of the one-loop effective potential 
at any temperature is gauge-dependent jiQll3^lia[l3l[ll[lll[l5l[l6l[171[l8l[ig 



so that the ratio on the left-hand side of Eq.(15) is not a well-defined physical quantity. 



Moreover, the procedure outline above to determine Tc perturbatively also introduces a 



spurious gauge-dependence, as one may observe by expressing the right side of Eq.(ll) 
in an arbitrary gauge [15] : 

PrT-'-^ . (16, 

where the additional terms are ^-dependent contributions associated with the one-loop 
corrections to the Higgs quartic self-coupling. As indicated above, the conventional 
analyses have been performed in Landau gauge, even though a small change in the 
choice of gauge parameter can significantly alter the ratio (pc/Tc. 

Obtaining a gauge-invariant baryon number preservation criterion (BNPC) requires 
several modifications of the naive perturbative treatment described above: 

1. Determine (or T„) in a gauge-invariant manner by following the evolution of 
Ves{4>,T) and consistently implementing the so-called Nielsen identities [HI 1^ . 
One procedure for doing so entails carrying out an /i-expansion of Ves{4>,T), as 
outlined in Refs. [121 HQ [12]. A generalization of this procedure also allows one to 
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approximate the full daisy resummation in a gauge-invariant way and to reproduce 
the trends for Tc found in nonperturbative calculations ^9] . 

2. Perform a gauge-invariant computation of the energy of the sphaleron configuration, 
Esph- In perturbation theory, it is possible to do so in the broken phase by working 
with the high-temperature effective theory in which zero-temperature masses are 
replaced by their (gauge-invariant) Debye masses. The energy £^sph then depends 
on a gauge-invariant scale v{T) that is not the same as 4>{T), and the fluctuation 
determinant k that characterizes the leading quadratic corrections to the sphaleron 
action [5n[52]. 

3. Compute the baryon density at the end of the EWPT, corresponding to a time 
delay of A^ew after its onset, and compare with the initial density resulting from 
the CPV transport dynamics described in Section [3] below. The resulting ratio is 
called the "washout factor" 

. (17) 

For the baryon asymmetry created by EWBG to be preserved, the washout factor 
S must not be too small. 

Rewriting the washout factor in terms of X, defined according to 

5 > e"^ , (18) 

the quantitative BNPC is [19]: 

4vr5 ^(^c)_^j^^^_j^^_j /AtEw\ (^g^ 



Here B parameterizes the relationship between the scale v{T) and the sphaleron energy 
[El [53] 

E^m^B^-^m\, (20) 



where S is a constant of order unity that depends on the mass of the Higgs boson, 
f (0) = 174 GeV is the value of the Higgs VEV at T = 0, and ayj is the weak 



coupling. The other quantities in Eq.(19) include: the Hubble time tn] a quantity Q 
characterizing the contribution of sphaleron zero modes; a function J-" that characterizes 
the dependence of the unstable mode of the sphaleron on v{T)-^ and a factor k accounting 



for fluctuations that are not zero modes. The appearance of the logarithms in Eq.(19) 
and the dependence on Atgw result from integrating the baryon number evolution 
equation [HI [M] 

UB , (21) 



dt 2 VT^ 

where Nf is the number of fermion families and T^^/VT^ oc exp{—Esph/T) is the 



sphaleron rate per unit volume inside the bubble. Qualitatively, the BNPC in Eq.(19) 
corresponds to the requirement that the sphaleron rate in the broken phase during the 
phase transition be much slower than the Hubble expansion rate. 
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Table 1. Maximum values of the Higgs boson mass, M^, for a first order EWPT in 
the SM as obtained from lattice studies. 

Lattice Authors (GeV) 

4D Isotropic [H] 80 ± 7 

4D Anisotropic [69] 72.4 ±1.7 

3D Isotropic [67] 72.3 ± 0.7 

3D Isotropic [65] 72.4 ± 0.9 



The conventionally-employed condition of Eq. ( 15 ) results from replacing the gauge- 



invariant ratio v{Tc)/Tc by the gauge-dependent one (f)c/Tc and making specific choices 



for the parameters appearing in Eq.(19). In particular, it has been assumed that one 
may take X = 10, corresponding to allowing the initial baryon asymmetry to be five 
orders of magnitude larger than what is presently observed - an assumption that is 
questionable in light of recent studies of the CPV transport dynamics discussed in 
Section [3j Additional significant uncertainties are associated with the value of the 
fluctuation determinant k and the duration of the transition At-^w- In short, even if one 
employs an appropriately gauge-invariant procedure to determine the degree of baryon 
number preservation, considerable uncertainty remains as to the precise requirement. 

Nearly all phenomenological studies carried out over the past decade or so have 
neglected these issues. Even if one places some trust in the use of perturbation theory 
to analyze EWPT dynamics relevant to EWBG, it should be clear that considerable 
work remains to be carried out in order to obtain a robust statement about the presence 
or absence of a sufficiently strong first order transition in a given BSM scenario. It may 
be, for example, that a given BSM scenario significantly modifies the dependence of Egph 
on the gauge-invariant scale v(Tc), the dependence of v on itself, the duration of the 
transition, or the fiuctuation determinant. From our view, then, conclusions that have 
been drawn as to the viability of EWBG based on existing perturbative analyses of the 
scalar field dynamics should be viewed as provisional at best and, in an ideal situation, 
revisited in light of these open theoretical issues. 

2.2. Non-perturbative studies 

Dedicated non-perturbative numerical analyses of the EWPT have been carried out for 
the SM, the Minimal Supersymmetric Standard Model (MSSM), and the two Higgs 

doublet model [SI [Ml [Ml EZl [Ml ESI EHl ED E21 E31 IMl ESI ^ 

Among the properties studied that are particularly relevant to EWBG are the critical 
temperature, T^; the weak sphaleron rate, Pgph (the focus of Refs. [571 [SS [ZQ]), and 
the character of the transition (first or second order, cross over, etc.). As we discuss in 



Section [4^ the latent heat L is also critical for the amplitude of gravity wave production, 
and it has been studied in, for example, Ref. [TU] . 

The process of bubble nucleation that is the starting point for EWBG requires 
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Figure 4. Dependence of the critical temperature Tc in the MSSM on the 
SUSY-breaking RH stop mass parameter mu taken from Ref. |66j (top panel) and 
Ref. [72] (bottom panel) . Note that the contribution to the tji mass-squared from 
this parameter goes as — = mfj . The phase labeled "CB" or "broken U" 
denotes a phase in which the stop field acquires a non-zero vacuum expectation value, 
corresponding to a color and charge-breaking vacuum. Top panel reprinted from M. 
Laine and K. Rummukainen, |Nucl. Phys. B 535, 423 (1998)| with permission from 
Elsevier. Bottom panel reprinted with permission from F. Csikor, Z. Fodor, P. Hegedus, 
A. Jakovac, S. D. Katz and A. Piroth, Phys. Rev. Lett. 85, 932 (2000)^ copyright 
2000 by the American Physical Society. 



a first order transition. In the SM, the presence of a first-order transition is found 
to require a sufficiently hght Higgs boson. Consequently, considerable attention has 
been paid to the value of the maximum value of the Higgs mass for which a first- 
order transition exists. Representative results from lattice studies are given in Table [T] 
The results obtained with the three-dimensional (3D) lattices require first carrying out 
the procedure of dimensional reduction to a the 3D effective theory (see, e.g., Ref. [61] ) 
before studying the phase transition properties of the latter using Monte Carlo methods. 
For either 3D or 4D studies, a criterion must be identified for determining the character 
of the phase transition. Among those employed are the susceptibility associated with 
the scalar field x ((0^ — (0^))^) correlation lengths. The scaling behavior of 
X with lattice volume can be used to determine whether the transition is first order, 
second order, or cross over. For Mh » ^ 75 GeV, as implied by collider searches 
for the Higgs, the transition appears to be a cross over transition. 
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A first-order transition can be accomplished by introducing additional scalar fields 
that couple to the Higgs sector. As discussed below, a particularly important example 
is the MSSM where the scalar top quarks couple strongly to the Higgs fields. Lattice 
studies imply that for sufficiently light right-handed (RH) stops, in, the transition can 
be not only first order but strongly so, and that decreases with the RH stop mass. The 
lowering of Tc is important for the protection of the baryon asymmetry from washout 
by EW sphalerons inside the expanding bubbles, since at the onset of the transition 
Tsph ~ Gxp{—Esph/Tc). Lowering Tc reduces Fgph by both reducing the magnitude of 
the denominator in the exponential and increasing the sphaleron energy -Egph- Results 
for the dependence of on the SUSY-breaking RH stop mass parameter obtained from 
Refs. [66l [72] are shown in Fig. |4] . Note that the stop mass-squared parameter must 
be negative to lower Tc. Increasing its magnitude in this parameter space region thus 
lowers the overall mass of the tn. 



2.3. Extending the SM scalar sector 

The need for a strongly first-order EWPT is one of the two reasons why EWBG does 
not work in the SM. Indeed, the failure of EWBG in the SM (or baryogenesis of any 
form) is a motivation for new physics near the electroweak scale. Such new physics is 
also needed to stabilize the electroweak scale itself, and can also account for the dark 
matter. While a broad range of extensions of the SM have been proposed to strengthen 
the EWPT to allow for EWBG, most of them fall into two groups. In the first group - 
exemplified by the MSSM - new scalars couple to the Higgs field and enhance the cubic 
term in the effective potential by running in loops. The second group consists of scalar 
fields coupling to the Higgs that develop non-trivial dynamics in the early Universe that 
infiuences the effective Higgs potential directly. In both cases, new hght scalars with 
significant couplings to the Higgs are needed, and these may lead to observable effects 
at coUiders. 

In the first class of new physics, in which a new scalar X modifies the Higgs 
potential through its loop effects, the relevant interactions can usually be written in 
the form [381 UH 

-£d M|.|XP + — |X|^ + Q|XniJ|2 . (22) 
6 

The third term is a Higgs portal coupling that cannot be forbidden by symmetries. 
Assuming that X does not develop an expectation value, the mass of the physical 
complex scalar is 

m\ = Ml + ^4>^. (23) 

Applying this to the effective potential yields 

AVe//(0,T) D [Yix{T) + Ml + Qct>'/2f" , (24) 

where Hx(T) is the thermal mass of X and ux is the number of degrees of freedom. If 
(50^/2 is much larger than the other terms for ~ 0^; this correction gives a strong 
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enhancement of the cubic operator that drives a first-order phase transition. If X is 
charged under SU{3)c, the contribution to the cubic is further enhanced at two-loop 
order by corrections involving virtual gluons [75l [76]. The net result is that a strongly 
first-order EWPT can obtain for Q > 1 and Mjr < if X is a SU{3)c triplet, but much 
larger Q values are needed if X is a gauge singlet [721 EI! • 

New scalars coupling in this way to the Higgs field have been studied extensively in 
the context of supersymmetry, and especially within the MSSM. For early perturbative 
studies, see Refs. [THl [79l [HD]. In this case the X scalar corresponds to the lightest 
scalar top quark (stop). It is a SU{3)c triplet with electric charge q = 2/3, and it 
must be mostly right-handed to be both light and consistent with precision electroweak 
data [38] . The other (mostly left-handed) stop must be considerably heavier to push up 
the mass of the Higgs boson [3H]. The couplings in Eq.(22) are fixed by supersymmetry 
to be Q — Ut and K ~ AttOs (up to small corrections from electroweak effects and 
left-right stop mixing), while the M]^ mass term corresponds to the right-handed stop 
soft mass-squared. 

A first-order phase transition strong enough for EWBG is found to be marginally 
consistent with the MSSM [3HI IMl EI]- As discussed above, in order to quench the 
sphaleron rate by lowering Tc, a negatively- valued Mx mass term is needed to cancel 
against the thermal mass appearing in Eq.(12) [3H]- This has the effect of destabilizing 
the origin of the scalar field space in the X direction, making possible an expectation 
value for the X field. Such a VEV would be disastrous since it would imply the 
spontaneous breakdown of electric charge and color. However, a detailed computation of 
the thermal evolution of the effective potential for and X finds that slightly negative 
values of M"^, but still large enough for EWBG, are phenomenologically acceptable. In 
this scenario, thermal effects can stabilize the X direction more strongly than the (j) 
direction, and the system can evolve first to a local minimum with 7^ and X = 0. 
This local minimum might not be as deep as the charge-breaking minimum with X 7^ 0, 
but it can be sufficiently long-lived to describe the Universe we observe. A negative 
also implies a very light stop mass, which forces the other heavier stop to be extremely 
heavy to drive up the mass of the lightest Higgs boson. 

A second way to make the electroweak phase transition more strongly first-order 
is to couple the Higgs to a new scalar that develops a VEV near the electroweak 
scale [82], |83l [H] . A simple example of this is 

-CD mlN^ + AnN^ + XnN^ (25) 
+ iAHN + CHN'')\H\^ + ... 

These interactions can allow both H and N to develop expectation values, resulting in a 
mixing between the physical singlet and SU{2)l scalars in the theory. When the singlet 
and SU (2) i mass parameters are similar, it is convenient to track the evolution of the 
VEVs in polar coordinates [82t 183] : 

{H^) = (f cos a, {N) = !fsma. (26) 
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Applying this parametrization to Eq.(25), one obtains cubic terms in the tree- level 
potential for ip that can lead to a strongly first-order electroweak phase transition. The 
singlet can also strengthen the phase transition by contributing to the loop-induced 
cubic term in the effective potential or by reducing the effective Higgs quartic coupling 
near the critical temperature [H3]. Similar effects have been found in gauge extensions 
of the SM [85] as well as theories with two or more SU {2)l doublets [731 EHl EH 188] . 

When the characteristic mass scale of the singlet sector is significantly larger than 
the SU (2) L part, the singlet can be integrated out of the theory. This produces effective 
Higgs interactions of the form 

- CD -fi^\H\^ + \\H\^ + ^\H\^ + ... (27) 

where A^v characterises the mass scale of the singlet sector. For values not too large, 
Atv < 1000 GeV, the new sextic term can also drive a strongly first-order electroweak 
phase transition [891 [9Ql [91] . 

Both classes of mechanisms to enhance the strength of the electroweak phase 
transition require new physics below the TeV scale. This is precisely the energy regime 
that is currently being probed directly by high energy collider experiments and indirectly 
by lower-energy precision probes. Some of the resulting experimental signals will be 
discussed in Section HI 



3. Creating CP asymmetries and baryons 

Baryons are created in EWBG through a combination of scattering, diffusion, and 
transfer reactions in the vicinity of the expanding bubble walls during the EWPT. 
Sphaleron transitions are rapid in the symmetric phase outside the bubbles, and they 
provide the requisite source of baryon number violation. However, sphalerons alone are 
not enough; a net CP asymmetry is also needed to bias the sphalerons to to produce 
more baryons than antibaryons. 

The CP asymmetry relevant for EWBG is n^,, the excess number density of left- 
handed fermions over their antiparticles. Such an asymmetry can be created during the 
EWPT together with an equal opposite asymmetry in riR by CP-violating interactions 
in the bubble wall. However, since sphalerons correspond to transitions between distinct 
SU (2) L vacua, only til directly affects baryon creation. The corresponding equation for 
the baryon density is 

d^fs = W^Kt, x)nB{x) + fcl2)(T, xWix)] (28) 

where = {nB,jB) is the baryon number current density, x is the coordinate orthogonal 
to the bubble wall, and k^I{T,x) are the weak sphaleron rate constants that account 
for the change in rate outside and inside the bubbles. Outside one has [M] 

A;(:)(r,x) =7^x|^ (29) 

out V I-^ 



X 
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with Tyjs/V giving the weak sphaleron rate per unit volume for Nf fermion famihes 

r 



n 



out 



15 



(30) 



and K 



20 ^54J. Deep inside the bubble these expressions become 



2 VT^ 







(31) 



(32) 



The computation of T^s/V inside the bubble has been discussed in Section [2] Note that 



Eq. (21) used in determining the washout factor is an approximation to Eq. (28) that 



neglects jB inside the bubble. It is clear from Eq.(28) that without the chiral asymmetry 
Hi, no net baryon charge will be created. 

The production of rii arises from interactions of particles with the bubble wall. 
The corresponding interaction rates are typically much faster than k^s{T,x) so that 



in practice one may decouple Eq.(28) from the system of equations that govern rii 



production. The latter encode a competition between several distinct dynamics: 

(a) C- and CP-violating interactions with the bubble wall that lead to the generation 
of particle number asymmetries for one or more species 

(b) particle number-changing reactions that tend to drive the plasma toward chemical 
equilibrium; 

(c) flavor oscillations that result from off-diagonal mass-matrix elements 

(d) scattering and creation-annihilation reactions that cause diffusion of the particle 
asymmetries ahead of the bubble wall and that push the system toward kinetic 
equilibrium in this asymptotic region. 

In the earliest work on this problem, these equations were reasonably assumed to 
have the form motivated by kinetic theory [22] 

d,Jk = - E r'V^ + S^""" , (33) 



al 



where = {rik^jk) denotes the current density for a given species of particle "fc", V 
is the rate for a given reaction "a " that affects the chemical potentials of particle 
species , and S'^^^ is a CP- violating source term. Diffusion arises in this formalism 
by making the diffusion ansatz, which relates the three-current to the gradient of the 
number density: 



where Dk{x,T) is the diffusion constant for the particle species. 



(34) 



The reaction rate terms in Eq.(33) couple different particle currents to one 



another. For example, the top Yukawa interaction involving the third-generation SM 
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quark doublet Q, the right-handed top quark T, and the Higgs doublet H leads to a 
term of the form 

- Ty {flQ - /iT - /iff) (35) 

in the transport equation for Jq that couples to and j^. CP-violation in this example 
is embodied by the source Sq^^ , that tends to generate a non-zero number density uq. 



This number density would subsequently be diluted by the reaction term of Eq. (35) 
that transfers some of it to the T and H densities. 



Under the assumptions leading to Eq.(33), a robust determination of n^^ entails 



identifying all of the relevant particle species involved in the plasma dynamics, 
computing the relevant CP-violating sources and particle number changing reaction 
rates, and solving the resulting set of coupled equations. Significant progress has been 
achieved over the years on all aspects of the problem. We first review studies of the 
particle number changing reactions. Next, we discuss recent advances in deriving reliable 
CP-violating sources. Throughout the discussion, we will use the MSSM as our primary 
example since this is the theory in which these issues have been studied in the greatest 
detail. 

Before proceeding, we make a few preliminary comments regarding timescales. An 
important consideration in determining the relative importance of a given process in the 
production of n/, is the associated time scale compared to the time scale for diffusion 
ahead of the bubble wall. As outlined in Refs. [H31 MM ES], in order for a given particle 
density to affect the dynamics of n^, generation in the unbroken phase, that density must 
diffuse ahead of the advancing wall. The time it resides in the unbroken phase is the time 
it takes for the advancing wall to "catch up" with the diffusing density. For a given time 
t, the diffusion length is d^ig = \/Dt, where D is an effective diffusion constant formed 
from an appropriate combination of individual diffusion constants Dk- The distance 
traversed by the wall is (i„aii = Vwt. Equating the two leads to the diffusion time 
scale: r^is = D/v"^. Taking representative estimates for the effective diffusion constant 
D ~ 50/T and wall velocity — 0.05 leads to T^is ~ 10^/T. Any process X having 
Tx « Tdig must be included in the analysis of dynamics leading to production, 
whereas for tx » T^is the process X effectively decouples. In particular, the rate for 
electroweak sphaleron transitions implies tew » ^difr so that one may first solve the 
dynamics leading to rii and then use the latter as an input for the generation of ub via 



Eq. (28). 



Note that the wall velocity plays a significant role through its impact on r^ig. As 
Vw increases, other processes have less time to equilibrate before the relevant particle 
species are captured by the bubble. This effect is particularly important when is 
sourced by CPV interactions involving BSM degrees of freedom. In order for any charge 
asymmetry involving the latter to be transferred into a non- vanishing ni , the particle 
number-changing reactions that facilitate this transfer must take place more quickly 
than the time it takes for the wall to overtake the diffusion process. 
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Figure 5. Particle number profiles with respect to the bubble wall {z — 0) obtained 
in the MSSM |95j . Pink, green, and red curves give, respectively, number densities 
of third generation left-handed quarks, third generation left handed leptons, and total 
TT-L- Bubble wall interior corresponds to z > 0. Thin curves represent results of an 
analytic approximation valid sufficiently far in front of the wall. Reprinted from D. J. 



H. Chung, B. Garbrecht, M. .J. Ramsey-Musolf and S. Tulin, |JHEP 0912, 067 (2009) 
with permission from JHEP. 



3.1. Particle number changing reactions 

Tracking the full set of coupled particle number densities can be very complicated. For 
instance, in the MSSM with CP-violation in the Higgsino-gaugino sector, the resulting 
set of coupled equations can involve up to 30 or more components. Fortunately, in many 
cases the scope of the problem can be significantly reduced. 

The most important simplification to the coupled equations is to relate algebraically 
the number densities of different particle species linked by reactions that are very 
fast compared to the other time scales of the problem (such as the diffusion time or 



the passage time of the bubble wall). In our top Yukawa example in Eq.(35), the 
corresponding reaction is typically very fast with ry^ << r^is, implying that this 
process achieves chemical (as opposed to kinetic) equilibrium well before the Q, H, 
and T densities are captured by the advancing wall. Consequently, it is an excellent 
approximation to set 

fXQ - fXT - fJ'H = ■ (36) 

This reduces the number of coupled equations by one. 

A frequent simplification arises in supersymmetric theories such as the MSSM as 
a result of "superequilibration" , wherein fip = fip with P and P denoting a Standard 
Model particle and its superpartner, respectively. Superequilibrium may arise either 
when a given supergauge interaction is fast compared to the timescale for diffusion 
ahead of the advancing bubble wall, or a chain of Yukawa reactions effectively yields 
/ip = /ip. Exceptions may occur, however, when either the gauginos become heavy, 
thereby suppressing the corresponding supergauge reaction rates, or a chain of Yukawa 
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reactions is kinematically blocked. 

Recently, the authors of Ref. |95] carried out an analysis of the particle number 
changing reactions for the MSSM under the assumptions leading to the general form 
of the transport equations in Eq.(33), obtaining numerical solutions to the full set of 
coupled transport equations. Superequlibrium is found to obtain in significant regions 
of the parameter space. Typical results for various particle number density profiles 
obtained in this work are given in Fig. [5] Baryon number generation is governed by the 
size of the diffusion tail for ul (red curve) ahead of the advancing wall < 0). 

For scenarios involving more than a single Higgs doublet H, such as the MSSM, 
an additional situation involving Yukawa interactions may arise that can substantially 
alter the generation of n^. In the Standard Model, the hierarchy of fermion Yukawa 
couplings implies that only ry^ is much shorter than rjifj and that one need consider 
only the Q—T—H interactions indicated above. In a two Higgs doublet model (2HDM) 
scenario, however, the bottom quark (squark) and tau (stau) Yukawa couplings may 
be enhanced with respect to their SM values by tan /3, implying that the Fy involving 
these particles grows as tan^/3. For tan/3 > 5 (20), the bottom (tau) Yukawa reactions 
effectively equilibrate on a timescale shorter than the diffusion time scale, implying that 
one cannot neglect their effect on the coupled set of transport equations [93] . 

In the MSSM case, the net effect of the enhanced bottom Yukawa coupling can be 
to significantly suppress the final baryon asymmetry ub/ s when the masses of the right- 
handed top and bottom squarks are nearly degenerate, or even lead to a sign change in 
the baryon asymmetry for m^^ < mi^. For sufficiently large and light tau sleptons, a 
significant til (and, thus, ub/s) asymmetry can nevertheless be generated, as an initial 
Higgs-Higgsino asymmetry is transferred into the left-handed tau lepton sector through 
the tau Yukawa reactions [931 El]- These dynamics are illustrated in Fig. [6| where one 
sees the vanishing of the baryon asymmetry for degenerate right-handed tops (ti) and 
sbottoms (6i). 



3.2. CP-violating sources 

As the foregoing discussion makes clear, analyzing the detailed nature of the particle 
number changing interactions and their parameter dependence is essential for predicting 
the baryon asymmetry, even if the CP-violating sources have maximal strength. At 
the same time, deriving reliable CP-violating sources remains an on-going theoretical 
challenge that is part of a broader quest for a systematic treatment of the EWBG 
transport dynamics. 

The earliest papers derived S'^^^ using computations of quantum mechanical 
reflection and transmission from a barrier (the bubble wall). The idea was originally 
applied to the generation of a lepton number asymmetry during the EWPT by Cohen, 
Kaplan, and Nelson [961 EZ] and subsequently applied directly to the baryon sector 
by Farrar and Shaposhnikov [HHl EH [100] and others (see, e.g. Refs. [1011 11U2] ). The 
resulting baryon asymmetry within the Standard Model, computed using this approach. 
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Figure 6. Baryon asymmetry as a function of right-handed bottom squark mass in 
the MSSM, scaled to the value obtained from the cosmic microwave background IHS] 
The blue (green) curves correspond to right handed stau mass of 1 TeV (90 GeV), 
while black dotted line is obtained neglecting the bottom and tau Yukawa couplings. 
Solid curves give results of full numerical solution to transport equations, while dashed 
line corresponds to an analytic approximation. A right-handed stop mass of 200 GeV 
has been assumed. Reprinted with permission from D. J. H. Chung, B. Garbrecht, M. 
J. Ramsey-Musolf and S. Tulin, Phys. Rev. Lett. 102, 061301 (2009), , copyright 2009 
by the American Physical Society. 

was generally considered to be too small to account for the observed value [221 123], so 
attention turned to scenarios involving physics beyond the Standard Model (BSM). A 
significant advance appeared in the incorporation of diffusion [92], with the resulting 
framework remaining the state-of-the-art for the better part of a decade when applied 
to BSM models (see, e.g. [TSl |79l |103l [lOSl [M [TUT]). 

It was subsequently realized that the essentially Markovian nature of this framework 
omitted potentially important "memory effects" that could lead to further enhancements 
of ul |1U8[ I1U9[ IllUj . Using the VEV insertion approximation depicted in Fig. [7| 
where the Higgs VEV-dependent off-diagonal entries of the mass matrix are treated 
as perturbative interactions, Ref. jllOj computed S""^^ in the MSSM for Higgsino and 
top squark transport near the bubble wall. As shown in the figure, the inputs for 
these sources involve two distinct particle species, either the Higgsinos and electroweak 
gauginos in the fermionic source case or the left- and right-handed top squarks in the 
scalar case. When the masses of the two species become nearly degenerate, S^^^ 
becomes resonantly enhanced. The presence of this feature is phenomenologically 
important, as an enhanced source will generate a particle asymmetry more effectively 
for a given CP-violating phase. Given the stringent limits on CP-violating phases 
implied by the non-observation of the electric dipole moments (EDMs) (see below), 
the resonant source enhancement is essential for viable EWBG in many instances. The 
same resonant enhancement was also found to apply to the CP-conserving parts of 
Fig. [7| which increases the relaxation rate of the non- vanishing densities sourced by the 
CP-violating parts, reducing the strength of the resonant enhancement of the baryon 
asymmetry but not completely mitigating it [lllj. 
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Figure 7. Contributions to CP-violating sources in the VEV-insertion 

approximation flllj . Here, v{x) denotes appropriate combination of scalar background 
fields (VEVs) at spacetime position a;^ interacting with (a) fermions /i,2 (b) scalar to 
quarks Il^r- Reprinted with permission from C. Lee, V. Cirigliano and M. J. Ramsey- 



Musolf, |Phys. Rev. D 71, 075010 (2005) , copyright 2005 by the American Physical 
Society. 



A rigorous treatment of the spacetime evolution of particle densities that 
incorporates the memory effects is properly made using the Schwinger-Keldysh Closed 
Time Path (CTP) formulation of non-equilibrium quantum field theory [TT2l 11131 1114[ 
1115] (for a pedagogical discussion in the context of EWBG, see Ref. |lllj ). In contrast to 
the "in-out" matrix elements relevant to scattering processes, analysis of the spacetime 
evolution of densities in the plasma requires study of the "in-in" matrix elements, 
appropriately averaged over a thermal ensemble. The corresponding evolution involves 
four Greens functions 

where for a complex scalar field (pix) 

G>{x,y) = {<P{x)ct>\y))T (38) 
G<(x,y) = (0t(y)0(x))r (39) 

while G\x,y) = 9{xo - yo)G>{x,y) + 9{yo - Xo)G<{x,y) and G\x,y) = 9{xo - 
yo)G^{x, y)+9{yo—xo)G^{x, y). Here, the subscript T indicates an appropriate ensemble 
average. The fully interacting Greens functions G{x, y) satisfy a pair of Schwinger-Dyson 
equations 

G{x, y) = G(x, y)o + I d^zG{x, z)otl{z, y)G{y, z) (40) 



G{x,y) = G{x,y\ + j d^zGix, z)Il{z,y)G{y, z)o (41) 

where the "0" subscript indicates the free Greens function and 11 is a matrix of self- 
energy functions corresponding to Eq.([37)). Acting with the free equation of motion 
operator ( e.g., + w? for a scalar field) on Eqs. ( 40||4T ), taking the difference of the 



two resulting equations, and considering the "<" component in the limit x ^ y yields 
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the transport equation 

dh / U>{X,z)G<iz,X) -G>iX,z)n<{z,X) 

J — oo 

+ G<{X,z)U>{z,X) -U<{X,z)G>{z,X)] , (42) 

where the effect of all interactions are contained in the convolution of the Greens 
functions and self-energies on the right hand side. In the VEV insertion approximation, 
CP-violating and CP-conserving interactions of the scalar field with spacetime- 
varying Higgs background field (the bubble wall) are contained the self-energies. 
Particle number changing interactions, as well as those associated with scattering and 
creation/annihilation, also live in the IT. Note that the source for the divergence of the 
particle number current contains an integral over the history of the system, leading to 
the presence of memory effects as emphasized in Ref. |110j . 

Subsequent work has attempted to refine this formulation, particularly within 
the MSSM. The authors of Refs. |116l 1117] endeavored to go beyond the VEV 
insertion approximation by noting that, when resummed to all orders in the VEVs, the 
interactions of Fig. [7] lead to spacetime-dependent mass matrices for the supersymmetric 
particles. In the case of stops, for example, the corresponding term in the Lagrangian 
is 



where 



={tl t*^ \ ' r ^ (43) 




Mll{x) =Ml + ml{x) + At{x) (44) 
Mrr{x) = M| + ml{x) + At(x) (45) 
M1p^{x) = ytv{x) [At sm. (3 {x) — jj, cos (3 {x)] . (46) 

Here, Mq and M-p are, respectively, the spacetime-independent third generation 
left- and right-handed stop mass parameters; mt{x) = ytv{x)/\/2 is the spacetime 
dependent top mass, with top Yukawa couphng yt and spacetime-dependent VEV 
v{x) = \Jv1^{x) + v1{x)] Ai(x) a function of /3(x) = ia.'nr^[vu{x)/vd{x)], v{x), and 
the weak mixing angle; and where /i and ytAt are supersymmetric Higgs mass and 
stop triscalar couplings, respectively. Note that in general Arg(/iy4t) ^ 0, so that 
the off-diagonal entry Ml^{x) effectively contains a spacetime-dependent CP- violating 
phase. At each point in spacetime, then, the propagating degrees of freedom - the mass 
eigenstates tj (j = 1,2)- are related to the weak interaction eigenstates ta (a = L, R) 
by a unitary transformation U{x) 

i, = U{x)]Ja . (47) 
The authors of Refs. |116[ 1117] computed the current densities for the ta (a = L, R) 



by expanding the mass-squared matrix M^^(x) in Eq.(43) about a given point z^ by a 



radius r, diagonalizing M^^{z) using Eq.(47), solving for the locally-free Greens functions 
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G^{z — r/2,z + r/2), and treating the first order term in the expansion of the mass- 
squared matrix, 

'1^ • («' 

as an interaction in the CTP Schwinger-Dyson equations, Eqs.(40 41). The resulting 
solution for the LH stop current j^^ then drives the creation of non-vanishing baryon 
number via the EW sphalerons, as per Eqs. ( 33p8 ). Similar methods were apphed to 
determine the chargino current j|±. 

The effect of including the fully re-summed VEVs and the first order interaction 
of Eq.(48) is to reduce the strength of the resonance. As an illustration for the 
chargino sources with M2 = 200 GeV, Ma = 150 GeV, and ATg{fj,At) = 7r/2, one 
finds that for ^ M2 the baryon asymmetry Yb is enhanced by about a factor of 
seven relative to its off-resonance magnitude at \fi\ = 100 GeV. This result represents 
about a factor of two reduction in Yg compared to the result obtained using the VEV 
insertion approximation in Ref. [78], and is considerably smaller still than obtained 
in Refs. |110[ 1111] . Nonetheless, it suggests that obtaining the observed value of Yb 
in the MSSM via chargino and neutralino sources would be viable even when a VEV 
resummation is performed. 

Parallel efforts on the resummation problem were carried out in Refs. 1120] , 

adopting the framework of the Kadanoff-Baym equations |121j that are deduced from 
the Schwinger-Dyson equations. (For a related effort using the Bogoliubov approach, 
see Ref. |122j ). Applying this analysis to the MSSM and performing a 2-fiavor analysis 
with approximate treatment of the off-diagonal densities the authors of Refs. |123t 1124] 
concluded that a proper treatment of flavor oscillations as well as a VEV resummation 
vastly suppresses Yb in the MSSM, rendering it unviable in light of current EDM limits 
on the CP- violating phases. 

To explain the arguments behind these conclusions, the subsequent response by 
the authors of Refs. |125[ 1126] . and the parallel work by the authors of Refs. |127t I128[ 
11291 1130] . we present here a specific example involving scalar particles, which avoids 
the complications associated with spin. Our example consists of a two-flavor scalar field 
theory. In the presence of a space-time dependent mass-squared matrix, the performing 
the flavor rotation of Eq.(47) leads to the Lagrangian in terms of local mass eigenstates 

= (01,02) 

^Cscalar = d^^^d^^ - ^^M''{x)<^ - ^tS'^a^^ (49) 

where M^(x) = diag[mf (x), m2(x)], 

= U\x)d^'U{x) , (50) 

and Cint contains interactions of $ with other particles in the plasma that lead to 
chemical and kinetic equilibrium as well as diffusion ahead of the bubble wall. 

The effects of CP-violation are encoded in E'^ that contains, in particular, a 
spacetime derivative of the effective CP-violating phase associated with the off-diagonal 
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Figure 8. Baryon asymmetry in the MSSM. Top panel: computed using VEV- 
insertion approximation and relaxation terms while including contributions from both 
charged and neutral Higgsinos and electroweak gauginos |lll) . Bottom panel: including 
only charged Higgsino and wino contributions but following the VEV resummation 
approach of Refs. |116l I117j . from Ref [131) . Top panel reprinted with permission 
from C. Lee, V. Cirighano and M. J. Ramsey-Musolf, P hys. Rev. D 7 1, 075010] 



(2005) , copyright 2005 by the American Physical Society. Bottom panel reprinted 
with permission from C. Balazs, M. S. Carena, A. Menon, D. E. Morrissey and C. E. 
M. Wagner, Phys. Rev. D 71, 075002 (200 5) , copyright 2005 by the American Physical 
Society. 



elements of the un-diagonalized matrix M'^{x)ab- Applying the equation of motion 



associated with all terms in Eq.(49) except Cint to the Schwinger- Dyson equations for 
the G leads to 



—I 



U{x,z)G{z,y) 



(51) 



with a similar equation holding for the operators acting on the y-coordinate. Taking the 
sum and difference of these two equations yields the constraint and kinetic equations, 
respectively. The former determines the dispersion relations for the propagating degrees 
of freedom, while the latter governs their dynamics. Extensive reviews of the related 
formalism can be found in Refs. |132t I133j . 

For the sake the present discussion, we concentrate only on the kinetic equation. 
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In doing so, it is convenient to consider the kinetic equation for the Wigner transform 
of the G^{x, y): 



G<{k,X) = j d^re'''-^ G<{x,y) (52) 

with r = X — y and X = {x + y)/2. The kinetic equation is 

2k-dxG-^{k,X) = e-^^|-z M\X) - 2ik ■ J:{X) + E\X),G<{k, X) 

+ C{k,X)^ , (53) 

where C{k,X) denotes the so-called "collision term" involving products of the G and 11 
Wigner transforms and the "diamond operator" is defined by 

I fdAdB dBdA\ 
oi^B) = - ——-——] . (54) 



2 XXt' k^ k„ 



Note that Eq. ( 53 ) and its partner constraint equation represent a set of coupled integral- 
differential equations in the space of Greens functions for the local mass eigenstates. For 
simplicity we have suppressed the corresponding indices that would appear, for example, 
in the term 



-2ik-J:iX),G<ik,X)]^^ = -2ik^ S^(X),,G<(A;,X) 



G<{k,X),,W{X),^ . (55) 



The different terms in Eq.(53) embody various key physical dynamics that govern 
the creation of a net particle density ahead of the bubble wall: (a) the [M^, G^] 
term gives rise to flavor oscillations, much in analogy with those observed or neutrino 
oscillations; (b) the commutators involving the field S'^ contain the effects of the 
spacetime-dependent CP-violating phases that are essential for the generation of a non- 
vanishing number density; (c) the collision term C includes the effects of scattering and 
particle creation/annihilation that gives rise to diffusion ahead of the bubble wall and 
thermalization within the plasma, thermal corrections to the masses and widths of the 
propagating degrees of freedom, and the particle number changing reactions discussed 



above. Note that Eq.(53) is exact to all orders in the spacetime- varying Higgs VEVs, 
is implied by the X^-dependence of the mass-squared matrix, the field S^, and C. The 
X'^-dependence of M^(X) also gives rise to the "semi-classical force" term that was 
first derived in this context for fermionic systems in Ref. |118j and that generates a 
CP-violating source term even for a single flavor situation. 

Obtaining exact solutions of the kinetic and constraint equations is a daunting task. 
Recent progress has been achieved in obtaining approximate solutions by expanding the 
various terms in powers of different scale ratios appropriate to the plasma dynamics 
(see Refs. |1111I125(I126] for a discussion). The scales include (a) the typical quasiparticle 
frequency or deBroglie wavelength Lint that is intrinsic to its free propagation. Both 
scales are set by the plasma temperature (assuming local thermalization): ~ L-^^^ ~ 
T. (b) The wall thickness or the timescale over which the particle masses vary 
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appreciably at a local point in space as the bubble expands; (c) The typical frequency 
of flavor oscillations, cUqsc given by the difference of the local eigenfrequencies — 
or the associated time and length scales, Tosc ~ Losc/c ~ w"^; the length or time 
scales associated with the collisional interactions, Tcoi ~ Lcoi/c. The latter set the rates 
for diffusion, kinetic and chemical equilibration, and thermal damping (widths) of the 
quasiparticles; and (d)the chemical potentials associated with various particle species 
fij. Typically one finds the following ratios of scales to be significantly less than one: 

Lint -t'int ,r-r\ 

^w = -j— Ccoi = — e^ = - . (56) 

The parameter is particularly important for electroweak baryogenesis, as it must 
be non-zero in order for the coupled set of kinetic equations to generate a non- vanishing 
particle density. Since typically one finds ~ (9(20/T), implying that < 0.5. 
In this regime, an expansion in powers of - sometimes referred to as the "gradient 
expansion" - is justified, and most recent work has focused on solutions obtained to first 
order in this small parametei[§] In general, one also finds that ecoi << 1 and << 1, 
leading to well-defined approximations by expanding in these quantities. In contrast, 
the ratio of "intrinsic" and oscillation length scales, eosc = -^int/-^osc depends strongly on 
the input parameters of the underlying Lagrangian. In the case of resonant baryogenesis, 
Uosc is relatively small, and one finds that the resonant enhancement occurs when the 
corresponding length scale Lose ~ L^ >> Lint. In this case, one may also expand in 
eosc (though some terms must be resummed to all orders in Wqsc in order to maintain 
consistency with the continuity equation). 

The work of Refs. |123l I124j yielded the first solutions of the kinetic equations for 
the MSSM under these approximations that take into account interplay of both the fully 
resummed spacetime varying VEVs with the off-diagonal densities. Applying the power 
counting in ew indicated above, the argued that the off-diagonal contributions to the 
RHS of Eq.(55) for the diagonal densities G^(A;,X)jj arise beyond leading non-trivial 
order, allowing one to neglect their effect as CPV sources for the diagonal terms. The 
corresponding numerical results then indicated that the value of Yb in the resonant 
regime is substantially smaller than obtained in earlier work. Given the limits on CP- 
violating phases from EDM searches (see below), one would then conclude that EWB in 
the MSSM cannot yield the observed baryon asymmetry, contrary to the implications 
of earlier work. 

Recently, however, the authors of Ref. |126j pointed out, using a simpler schematic 
two-flavor scalar field theory, that the approximations used in Refs. |123t I124j do not 
consistently implement the power counting in the e^. In particular, the solution to 
Eqs. (53) involve a integral of the terms in Eq.(55) over a distance scale of order 1/L^. 
Although the fields S'^ are nominally 0{eu)), this integral compensates for this nominal 
suppression. Thus, all contributions from the commutator in Eq.(55), including 
those involving off-diagonal terms k ■ Si2G'^i — Cfg^ ■ S21 that contribute to the evolution 

§ This validity of the gradient or expansion corresponds to the "thick waU" regime in earher work. 
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of the diagonal Greens functions Gfi must be kept to leading non-trivial order in the ej. 
Dropping these off-diagonal contributions, as in Refs. |123[ 1124] , removes the dominant 
CP-violating source for the diagonal Greens functions that characterize the diffusion of 
particle number ahead of the advancing wall and leads to the suppressed asymmetry as 
obtained in that work. On the other hand, the consistent solution that retains these 
terms displays the substantial resonant enhancement that was noted in earlier works 
obtained with the VEV-insertion approximation. 

It is important to emphasize that the most recent work on this problem has 
been completed using schematic two-flavor scalar field models, and not the MSSM for 
which CP-violating dynamics involving fermion fields (Higgsinos and gauginos) is the 
most viable mechanism for successful baryogenesis in this scenario. In addition, the 
authors of Ref. |13U] - building on earlier work [127^ 1128] 1129] and concentrating on the 
time-dependent, spatially homogeneous background field problem - have argued that 
contributions from so-called "fast modes" or "coherence shells" may have a significant 
impact on the magnitude of the particle asymmetries. In contrast to solutions of the 
constraint equations that satisfy particle-like dispersion relations (quasiparticle modes), 
the coherence shells correspond to regions of vanishing space or time components of the 
momentum and generate singular contributions to the spectral densities. A consistent 
treatment of these modes, which can mix with the quasiparticle modes, appears to 
require a resummation to all orders in ew The implications for the full spacetime- 
dependent problem remain to be analyzed. On the other hand, the progress achieved 
in Refs. |125] 1126] represents a significant advance, not only for the implementation of 
a consistent power counting, but also for the full inclusion of the collision term C{k,X) 
that governs the evolution of kinetic equilibrium ahead of the advancing wall. Previous 
work had largely implemented equilibrium as an ansatz, rather than obtaining it as a 
direct solution of the dynamics. 

Significant formal and phenomenological work remains to be completed. A full 
application to the dynamics of fermions - including resolution of questions involving 
the "semi-classical force" term, the coherence shells, and off-diagonal densities - as well as 
a more realistic set of thermalizing interactions in the collision term are among the open 
theoretical problems. Phenomenological application to realistic scenarios like the MSSM 
that involve a more extensive set of coupled integral-differential equations, including 
both the network of particle number changing reactions as in Refs. [95l 193] 1134] as 
well as the consistent treatment of flavor oscillations and spacetime varying background 
fields in the presence of the thermalizing interactions in the plasma, will entail significant 
effort. Nonetheless, given the recent theoretical progress and indications that resonant 
EWB remains a viable mechanism, this investment of effort is likely to have a significant 
impact on the field. 
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4. Testing electroweak baryogenesis 

Electroweak baryogenesis requires new particles and interactions to obtain a strongly 
first-order electroweak phase transition and to provide sufficient CP violation. These 
new particles cannot be much heavier than the electroweak scale, and they must couple 
significantly to the Higgs field. Together, these properties imply that such new particles 
will lead to observable effects in upcoming high-energy and high-precision experiments. 

The prospect for observing new particles directly at the LHC and indirectly through 
high-sensitivity, low-energy studies of CP- violating observables imply that EWBG is a 
generally testable and, therefore, falsifiable, baryogenesis scenario. In this respect, it 
contrasts with other scenarios that typically involve higher scales, such as standard 
thermal leptogenesis or Affieck-Dine baryogenesis. In what follows, we summarize some 
of the primary experimental tests of EWBG, focusing largely but not exclusively on 
the MSSM as an illustration. We consider tests that may be performed each of the 
three-frontiers in particle physics: the high-energy frontier; the intensity frontier; and 
the cosmological frontier. 



4.1. The Intensity Frontier: CP-Violation and EWBG 

In general, the most powerful probes of BSM CP-violation that are relevant for EWBG 
are searches for the permanent electric dipole moments (EDMs) of the electron, neutron, 
and neutral atoms. In all cases, only null results have been attained to date, implying 
stringent constraints on new sources of CP-violation. Limits also exist on the muon 
EDM, though it is considerably less constraining than the electron EDM limit, even 
after accounting for the me/rn^ suppression of the former with respect to the latter. 
At the most basic level, one expects the one-loop EDM of an elementary fermion / 
generated by new field(s) of mass M to go as 

.,^e(^)2.„, (57) 

where ak is either the fine structure constant or strong coupling (evaluated at the scale 
M) and is a CP-violating phase associated with the new interactions. For = aem 
Eq.(57) gives 

"'-^'"^ (m^) (— j ■ (58) 

The present limit on the EDM of the electron, \de\ < 10.5 x 10~^^ e cm |135] obtained 
from an experiment on the Yb-F molecule, then implies that 

|sin0| < (M/2 TeV)^ . (59) 

Similar constraints follow from the limits on the neutron |136] and ^^^Hg atomic |137] 
EDMs: 

|ci^(i99Hg)| <3.1 X 10-29 ecm , (61) 
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Figure 9. Two loop EDMs in the MSSM (see Ref. |141j and references therein). 
Reprinted with permission from Y. Li, S. Profumo and M. Ramsey-Musolf, |Phys.| 
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assuming that any contributions from the QCD O-term interaction are sufficiently small 
that no cancellations between this source of SM CP-violation and that arising from new 
interactions occurs. Contributions from CP-violation associated with the SM CKM 
matrix first arise at three (four) loop order for dn and dA (de), implying effects at well- 
below the 10~^° e cm level. The next generation of lepton, neutron, and neutral atom 
EDM searches aim to improve the level of sensitivity by up to two orders of magnitude in 
the short term, while efforts to reach even greater sensitivity with storage ring hadronic 
EDM searches are underway (for a recent summary of present plans, see e.g., Ref. |138] ). 

The constraints implied in Eq.(59) generically render EWBG unviable. For the 
new particles to be sufficiently abundant in the electroweak plasma at T ~ 100 GeV, 
their masses should be lighter than ~ 500 GeV, implying | sin0| < 0.01. In this case, 
the CPV sources in the transport equations discussed in Section [3] are suppressed and 
EWBG becomes untenable as a result (cf Figs. 6]|8). There exist, however, several paths 
to evading the one-loop EDM constraints. In the MSSM, the one-loop EDMs contain one 
scalar [e.g. squark or slepton) and one fermionic superpartner (gaugino or Higgsino). By 
making one or the other species sufficiently heavy, the one-loop EDMs can be evaded, 
thereby relaxing the constraints of Eq.(59) on the CPV phase. At the same time, 
the other superpartner species may remain relatively light, enabling its interactions to 
generate the CPV sources for electroweak transport dynamics. Large differences in the 
scalar and fermion mass spectra have been motivated on other grounds recently, as in the 
case of "split SUSY" models |139t 1140] that contain heavy ffist and second generation 
sfermions but relatively light Higgsinos and electroweak gauginos. The present generic 
LHC lower bounds on the masses of the gluinos and first and second generation squarks 
is at least consistent with this scenario. 

From the standpoint of CPV for EWBG, one must still consider EDMs generated 
at two-loop level, as in the "Barr-Zee" graphs of Fig. |9l Recently, the full set of 
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Figure 10. Curves of constant de and cZ„ for MSSM EWBG in the limit of heavy 
sfermions [143] . For each curve, is set to the value giving the correct baryon 
asymmetry. Reprinted from V. Cirigliano, Y. Li, S. Profumo and M. J. Ramsey- 
Musolf, [jHEP 1001, 002 (2010)| with permission from JHEP. 



such diagrams were computed by the authors of Ref. |141] and the corresponding 
imphcations for MSSM baryogenesis dehneated in Refs. |142[ I143j . Even with the 
two-loop suppression and the most optimistic CPV sources computed using the VEV- 
insertion approximation, only one CPV phase remains sufficiently unconstrained as to 
remain a potentially viable driver of MSSM EWBG: the relative phase of the bino soft 
mass parameter Mi, the soft parameter b and the supersymmetric Higgs mass parameter 
yu: 01 = Arg(Mi/i6*). Since 0i is associated with the presence of the bino degrees of 
freedom, it only enters the two-loop graphs involving the Xk exchange of a 

{W^, H^) pair, representing a small sub-class of the full two-loop diagrams. During the 
EWPT, the relevant CP-violating sources would be those involving the Higgsino-bino 
processes in Fig. [7} corresponding to "neutralino-driven" EWBG. A summary of the 
relation between baryon asymmetry and EDMs of the electron and neutron are given in 



Fig. 10 in the region of resonant electroweak gaugino-Higgsino EWBG. As indicated by 
the inner most contours, to conclusively test or rule out MSSM EWBG would require 
improvements in the sensitivity of EDM searches by roughly two orders of magnitude, 
roughly consistent with the goals of the next generation of experiments. 

Going beyond the MSSM, it is possible to introduce new CPV interactions in the 
scalar sector that could evade present EDM constraints but still generate the CPV 
sources as needed for EWBG [Ull HHl [US]. The authors of Ref. [HI] studied the 
CPV sources in the NMSSM and found that the presence of the additional gauge singlet 
superfield gives rise to a new CPV source that is second order in the expansion - 
associated with a "semi-classical force" term in kinetic theory - and that may contribute 
strongly away from the resonant regime. This source depends on the same CPV phases 
as in the MSSM, so one must contend with constraints from EDM searches. At the time 
this work was completed, a minimum first and second generation sfermion mass of 1 TeV 
was sufficient to evade the existing d^^n bounds. More recently, the authors of Ref. |146J 
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observed that a new phase associated with the gauge-singlet extension of the MSSM 
could successfully drive EWBG through both top and stop sources as well as Higgsino- 
bino interactions (for a related study, see Ref. |147j ). The dominant constraints on this 
phase are associated with two-loop contributions to the down quark "chromo-electric 
dipole moment" as it might generate a ^^^Hg atomic EDM. Nonetheless, NMSSM EWB 
remains viable even with these stringent constraints. An earlier and more extensive 
U(l)' extension was studied by the authors of Ref. [HS], who found a sufficiently large 
baryon asymmetry could be generated fromCPV sources associated with spontaneous 
CPV during the EWPT and whose CPV effects relaxed to sufficiently small values in 
the broken phase so as to evade the EDM bounds. 

In all of the foregoing scenarios, the CPV interactions are flavor diagonal. Recently, 
some attention has focused on the possibility that flavor non-diagonal CPV might 
provide the source for EWBG while evading the one-loop EDM bounds [HH 11481 1149] . 
As with SM CKM CPV, flavor non-diagonal CPV interactions that involve the second 
and third generation quarks would not contribute to EDMs until at least two-loop order. 
On the other hand, the associated CPV effects could appear in the B- and/or D-meson 
system, allowing for a potential probe of this possibility. To illustrate, we consider the 
schematic 2DHM of Ref. |149j . applied to the second and third generation fermions as 
needed for CPV in the 6 — )■ s transitions. The corresponding interaction is 

-y-Q^H,u^ + yf^Q\eHi)di, + h.c. , (62) 

where run over the second and third generations. For purposes of illustration, one 
may make several simplifying assumptions, including taking tan /3 = 1 at T = and 
setting ysh = = and retaining the one CPV phase that remains in the limit 
of vanishing yss and Xss (after field redefinitions): = Arg(Afes). In the VEV- 

insertion approximation the CPV sources are proportional to {Xbsl"^ sinOx^^- On the 
other hand, CPV in the Bg system arises from tree-level exchange of the Higgs field 
Hfjs = — cos (3 Hu + sin /3 if], leading to the mixing operator 

i^{bLSR){bLSR), with Abs^mjj^Jv 

with ~ m^j^Jv and Cbs |Abs|(l ± exp(i6'Ai,J| depending on the signs of the yij. 

The resulting curves of constant baryon asymmetry, shown in the (sin 9\^^, Xbs/ V^) 



plane are given in Fig. 11, along with constraints from the Tevatron and LHCb 
experiments. Note that in this restricted, schematic model, non-negligible contributions 
to the baryon asymmetry can occur. Future improvements in the sensitivity of LHCb 
studies could either uncover CPV in a region consistent with a significant contribution 
to the asymmetry or place stringent constraints on this possibility. On the other hand, 
EDM searches are relatively insensitive to the CPV phase 6x^^, as the operator in Eq.(63) 
contains no first generation quarks and involves only fiavor non-diagonal combinations 
of the opposite chirality quarks. Clearly there exists considerable room for additional 
studies of fiavor non-diagonal EWBG, including a more extensive study of the framework 
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sign(^s)fts 

Figure 11. Contours of constant baryon asymmetry in the "Electroweak 

beautygenesis" scenario of Ref. |149j . Figure reprinted with permission from T. Liu, 
M. Ramsey-Musolf and J. Shu, Phys. Rev. Lett. 108, 221301 (20T2)| copyright 2012 
by the American Physical Society. 

discussed here as well as those analyzed in Ref. [571 I148j . The extent to which this 
possibility could provide a viable EWBG alternative to flavor diagonal CPV in EWBG 
in light of increasingly stringent EDM limits remains to be seen. 

4-2. The High Energy Frontier: CP-Violation and Phase Transitions 

New particles related to EWBG can potentially be created at observable rates at high- 
energy colliders such as the LHC. This is especially true of new colored states that help 
to strengthen the EWPT, such as a light stop in the MSSM or more generally a new X 
scalar as discussed in Sec. [2j To have an adequate effect, these particles must be lighter 
than about mx ^ 200 GeV implying very large LHC (and Tevatron) production cross 
sections. Even so, such states be consistent with existing collider limits. 

One way to hide a light colored scalar is to have it decay to dijets. This can arise 
from a XqiQj coupling. Even though the dijets in this case would form an invariant- 
mass peak at the X mass, the backgrounds at the LHC (and the Tevatron) are so large 
that they swamp the signal in the low-mass region. Moreover, the cuts used in existing 
searches single or paired dijet resonances at both the LHC and the Tevatron do not 
extend or have limited sensitivity to lighter masses mx ^ 200 GeV |150[ I15H 1152] . 
Some sensitivity could be recovered with more generous cuts jl53j . or if the X decays 
involve heavy flavor |154j . 

A second way for a light colored state X to have escaped detection is for it to 
decay to a light quark and a long-lived neutral N fermion (that could be the dark 
matter). This is natural in the MSSM, where the light stop {X = ti) can decay to a 
charm quark and the lightest neutralino (A^ = Xi)- small X — N mass splittings, 
the decay products are very soft and difficult to detect using the standard searches for 
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multiple jets and missing energy, and dedicated Tevatron searches only limit splittings 
above about 35 GeV j l55[ 1156] . In Fig, 12 we illustrate the current limits on a light 
stop that decays via ii — )■ cXi based on data from the Tevatron and the LHC. Instead, 
some sensitivity can be recovered through monojet searches, requiring a single hard jet 
and missing energy |157l 11581 1159[ 1160] , as can also be seen in Fig. 12 , although a full 
test of the parameter space will be challenging |161j . A potential further probe may 
be provided by the decays of stoponium, a bound state consisting of a stop and its 
anti-particle [IS21 UMl MM MS\ ■ 

The new physics required for EWBG must also couple significantly to the SU{2)l- 
doublet Higgs field, and this can potentially induce observable changes in the production 
and decay properties of the Higgs boson. A very significant effect arises from colored 
X scalars that couple to the Higgs as in Eq.(22). Such scalars will contribute in loops 
to the amplitudes for Higgs boson production through gluon fusion and Higgs decay to 
diphotons. Large positive values of the coupling Q > are needed to induce a strongly 
first-order phase transition. In this case their contribution to gluon fusion interferes 
constructively with the dominant top quark loop in the SM, and destructively with the 
dominant loop for diphoton decay. The net result is a significant enhancement in 
the rate of gluon fusion that is closely related to the strength of the electroweak phase 
transition, and a more modest decrease in the branching fraction to diphotons [71] . 
Similar results are found for the light stop of the MSSM [166j. 

In Fig. 13 we show the enhancement of the gluon fusion rate from X relative to the 
SM (red dotted lines) as a function of Q and the mass parameter M|^, together with 
an estimate for where the phase transition is strong enough for EWBG (to the right 
of the thick black solid line). Gluon fusion is the dominant Higgs production mode at 
the LHC, and is the primary production channel for the 77, W^W~ and Z^Z^ decay 
searches at the LHC. Indeed, The enhancement of the gluon fusion rate implied by 
this mechanism of strengthening the electroweak phase transition is already strongly 
constrained by current LHC and Tevatron Higgs searches [7^ I167j . 
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Figure 13. Rates of gluon fusion Higgs boson production relative to the SM (red 
dotted lines) and contours of (pc/Tc (black solid lines) for one new color-triplet scalar for 
given values of the parameters —sgn (M|.) a/|M|.| and Q. The thick solid corresponds 
to 4>c/Tc — 0.9 and the adjacent solid lines show steps of A{(j)c/Tc) — 0.2. Figure from 
Ref. [75. 



New uncoloured X particles coupling to the Higgs boson as in Eq.(22) can also make 
the phase transition more strongly first-order, although their effect tends to be weaker 
than the coloured case. If such a state has a non-trivial electric charge, it will modify the 
Higgs branching fraction to — )■ 77 [TU EH]- As in the coloured case, the interference 
with the loop is destructive when the phase transition is made more strongly first 
order. The net effect, therefore, is to decrease the Higgs branching to diphotons while 
leaving the gluon fusion rate largely unchanged. Uncoloured X particles could also 
potentially be probed directly by their electroweak production channels, or their effects 
on the Higgs self-coupling |168] . 

The electroweak phase transition can also be made more strongly first-order if there 
are other fields that develop VEVs in the early Universe at about the same time as the 
Higgs. A simple example of this is the real singlet model presented in Section [2] [83] . 
In this case, there will be an additional fundamental scalar boson in the theory that 
will mix with the S'?7(2)i-doublet Higgs excitation. The resulting real scalar mass 
eigenstates will consist of a SM-like hi and a singlet-like h2- Decays of hi are frequently 
similar to the SM, but can be changed radically if hi — )■ h2h2 is kinematically allowed. 
The decays of the h2 state are typically inherited from hi, so the chain hi — h2h2 is 
likely to produce 46, 262r, and 4r final states, which can be distinctive but challenging 
to find a hadron colliders |169l 1170] . On the other hand, the singlet-like /12 state may 
have more exotic decay channels if there are other light states in the theory, such as a 
light singlet fermion |171l 1172] . More generically, one could expect a significant signal 
reduction in conventional SM Higgs search channels due to the effects of mixing and 
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hi —7- h2h2 decays if kinematically allowed [83], as well as the appearance of the second 
state /i2 |173j . 

New light particles are also needed to induce CP violation in the expanding bubble 
walls. In many cases, they carry non-trivial electroweak charges and couple to the 
varying-Higgs background. This is true of the MSSM, where the main source comes 
from light neutralinos and charginos |174] . Direct searches for such particles created 
via their electroweak production modes are underway at the LHC, and some relevant 
exclusions have been obtained |175] . Even so, the detailed signals depend very sensitively 
on the decay channels of the new states. 

4-3. The Cosmological Frontier: Gravity Waves and More 

The strongly first-order phase transition needed for EWBG can produce a cosmological 
signal in the form of gravity waves[T75l 1177] . As discussed above, the phase 
transition proceeds by the formation of bubbles of the electroweak broken phase 
within the surrounding symmetric-phase plasma. Gravitational radiation is created 
by the turbulent expansion of the walls |178] and their subsequent collisions as they 
coalesce |176[ I177j . The net effect of the many bubble collisions that occurred within 
the current Hubble radius is a uniform stochastic background of gravitational radiation 
with a characteristic spectrum. 

The spectrum and intensity of gravity waves created by a strongly first-order phase 
transition depend on three parameters characterizing the transition |177l 1178] : the latent 
heat a released by the phase transition at the nucleation temperature T„ relative to 
the background radiation energy, the characteristic rate of bubble nucleation /3, and the 
bubble wall velocity Vb- All three quantities can be computed from the finite-temperature 
effective potential discussed in Sec. |2j 

Estimates of the gravity wave signals produced by a strongly first-order electroweak 
phase transition suggest that it will be very difficult to detect in the foreseeable 
future |179t 118m 1181] . The signal is typically too low in frequency to be picked up 
by the LIGO experiment, but it may be visible at LISA if the transition is extremely 
strong. The prospects of discovery are considerably better at BBO, but even in this case 
the signal from the phase transition could be obscured by larger gravity waves signals 
due to astrophysical processes or inflation |180] . 

New physics related to EWBG can also play an important role of other aspects 
of cosmology. A specific example of this is dark matter (DM) in the MSSM. Here, 
the lightest neutralino can play a key role in generating the baryon asymmetry as 
described in Sec.[3| and can also make up the DM if it is the lightest superpartner (LSP). 
When the LSP is Bino-like, as motivated by the limits from EDM searches, the relic 
abundance tends to be too large. However, MSSM EWBG also requires a very light 
stop, and this state can reduce the Bino-like LSP abundance to an acceptable level by 
coannihilation |182] 1131 ]. A similar interplay between the new states contributing to 
EWBG and those making up the DM is also found in variety of other models. 
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5. Summary 

Electroweak baryogenesis remains a theoretically attractive and experimentally testable 
mechanism to generate the baryon asymmetry of the Universe. During the past decade, 
much of the theoretical attention has focused on leptogenesis - motivated in part by 
the discovery of neutrino oscillations and the relative simplicity of this scenario. More 
recently, however, there has been a resurgence of interest in EWBG due to its testability 
and the advent of new probes of physics at the terascale. Searches for new particles at 
the LHC could discover degrees of freedom that were thermally abundant during the 
electroweak phase transition and whose interactions could have engendered a strong first 
order transition while giving rise to particle asymmetries needed to seed the production 
of baryons. At the same time, a generation of new searches for permanent electric dipole 
moments will provide ever more powerful probes of possible flavor diagonal CP-violation 
in these new interactions that is an essential ingredient for successful EWBG. In parallel, 
new mechanisms involving flavor non-diagonal terascale CP-violation are under study, 
and these mechanisms may have signatures in experiments involving B-mesons at the 
LHC or super-B factories. In short, one may anticipate that either the ingredients for 
successful EWBG will be uncovered during the coming decade or that this scenario will 
be sufficiently constrained that more speculative, high-scale baryogenesis scenarios such 
as leptogenesis are left standing as the most viable alternatives. 

Achieving this confrontation of experiment with EWBG requires having in hand 
the most robust theoretical tools for computing the baryon asymmetry and a sufficiently 
broad phenomenological framework. In both respects, the past decade has witnessed 
considerable advances. Substantial effort has been devoted to deriving and solving 
the relevant set of transport equations that underlie the production of left-handed 
particle number asymmetries, while new attention has focused on achieving more 
theoretically well-defined computations of the phase transition properties. Nonetheless, 
there exists ample room for further improvements. In the case of the EWPT, 
future work could include development of additional methods for carrying out gauge- 
invariant perturbative computations, calculations of the fluctuation determinant in BSM 
scenarios, and investigation of the other sources of uncertainty that enter the BNPC of 



Eq.(19). New Monte Carlo studies in representative BSM theories would also provide 
important benchmarks for gauging the validity of perturbative computations and their 
phenomenological implications. Continued reflnements of the transport machinery 
would include application of the recent developments for scalar flelds to fermionic 
systems, updated analyses of the "semiclassical force" terms, and resolution of questions 
surrounding the coherence shells. From the standpoint of phenomenology, an interesting 
new direction entails the possible role of flavor non-diagonal CP- violation. Further study 
of signatures of BSM scalar sectors, including modifled Higgs production cross sections, 
branching ratios, and exotic flnal states - as well as scenarios in which interesting scalar 
sector extensions could evade discovery at the LHC - is also an obvious priority. 

More broadly speaking, explaining the origin of the visible matter of the Universe 
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continues to be one of the primary motivations for seeking what hes beyond the Standard 
ModeL While it is by no means certain that the explanation lies at the terascale, 
the time is ripe to address this possibility with vigor as part of the larger effort to 
determine what new symmetries and degrees of freedom - if any - are associated with 
the electroweak chapter of cosmic history. Whatever the outcome of the EWBG theory- 
experiment interface, one can expect this endeavor to yield important new insights into 
the fundamental laws of nature and their cosmological implications. 
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